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ABSTRACT 

We present echelle spectroscopic data for five eclipsing binary stars and two giant stars 
in the field of the open cluster NGC 2243. The average cluster velocity is determined to be 
+60.4±0.6 km s —1 . Four of the eclipsing binaries are very likely members of the cluster based 
on their observed radial velocities. The absolute parameters of cluster member NV CMa are 
determined by analysing photometric and radial velocity data. We obtain 1.089±0.010Mq and 
1.221±0.031i? Q for the primary, and 1.O69±O.O1OM and 1.178±O.O37R for the secondary. 
Both components of the binary are located on the Main Sequence, about 1 mag. below the 
turn-off point on the cluster color-magnitude diagram. Using model age-luminosity and age- 
radius relations we obtain 4.35±0.25 Gyr for the age of NV CMa. The derived age is, however, 
very sensitive to the adopted metallicity of the cluster. We demonstrate that a meaningful 
determination of the ages of objects like NV CMa based on evolutionary models is possible 
only if their metallicity is know with a relative accuracy of a few percent. The distance 
moduli calculated for the components of NV CMa agree closely with each other, and imply 
an apparent distance modulus of the cluster of (m — M)v = 13.25 ± 0.08. 

Stars: binaries: eclipsing, binaries - stars: individual: NV CMa - open clusters and asso- 
ciations: individual: NGC 2243 



1 Introduction 



The field of the intermediate-age open cluster NGC 2243 contains 5 known de- 
tached eclipsing binaries. Four of these are likely cluster members based on their 
photometric properties (Kaluzny et al. 2006; hereafter KKTS). The components 
of these systems should have the same age, metallicity and heliocentric distance 
and the determination of their absolute parameters can provide an interesting 
test of evolutionary models of low mass stars. Moreover, through the use of 
the surface brightness method, one may obtain a direct measure of the cluster 
distance. 

This paper is focused on the determination of the absolute properties of 
NV CMa, an eclipsing binary located in the central area of NGC 2243. In 
addition we use spectroscopic observations to check on the membership status 
of 4 other binaries in the cluster field. 



"This paper includes data obtained with the 6.5-meter Magellan Telescopes located at Las 
Campanas Observatory, Chile. 
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2 Spectroscopic Observations and Reductions 

Spectroscopic observations of NGC 2243 stars were carried out with the MIKE 
echelle spectrograph (Bernstein et al. 2003) on the Magellan II (Clay) telescope 
of the Las Campanas Observatory. The data were collected during observing 
runs in 2004 October and 2005 September. For this analysis we use data ob- 
tained with the blue channel of MIKE covering the range from 400 to 500 nm 
with a resolving power of A/AA«38,000. All of the observations were obtained 
with a 0.7 x 5.0 arcsec slit and with 2x2 pixel binning. At 4380 A the reso- 
lution was ^2.7 pixels at a scale of 0.043 A/pixel. The seeing ranged from 0.7 
to 1.1 arcsec. The spectra were first processed using a pipeline developed by 
Dan Kelson following the formalism of Kelson (2003, 2006) and then analysed 
further using standard tasks in the IRAF/Echelle package^. Each of the final 
individual spectra consisted of two 600 s exposures interlaced with an exposure 
of a thorium-argon lamp. We obtained 12 spectra of NV CMa. The average 
signal-to-noise ratios range from 25 at shorter wavelengths to 50 at longer wave- 
lengths. In addition to observations of NV CMa we also obtained single spectra 
of the eclipsing binaries V4=NS CMa, V5=NX CMa, V7 and V9 (names of vari- 
ables as in KKTS), and for the red giants H3110 and H4115 (Hawarden 1975) 
located in the central part of the cluster field. 

2.1 Spectroscopic Orbit of NV CMa 

Radial velocities of the components of NV CMa were measured by cross-correlation 
with the FXCOR task in IRAF, using observations of HD 33256 as a template. 
According to Nordstrom et al. (2004), HD 33256 has V rad = 10.1 ± 0.2 km s" 1 
and a projected rotational velocity Vsmi = 10 km s _1 . With B — V = 0A5 and 
[Fc/H] = —0.5, it provides a good match for the color index and metallicity of 
the binary. The template was observed with the same instrumental configura- 
tion as the variable. The correlation peaks were measured to have a FWHM of 
about 70 km s _1 . All spectra showed the peaks to be separated by more than 
1.4 FWHM, and the peaks were measured simultaneously with the FXCOR 
package. The correlation was measured only from the metal lines, excluding 
the H/3, H7, and H<5 hydrogen lines. Our velocity measurements for NV CMa 
are listed in Table 1 and are shown in Fig. 1. A Keplcrian orbit was fitted to 
the observations by fixing the period and epoch based on the precise ephemeris 
established by KKTS: 

Mini = HJD 2448663.70748(5) + 1.18851590(2) (1) 

We assumed a circular orbit based on the photometric data. This assumption is 
further supported by the relatively short orbital period of the variable compared 
to the cluster age of about 3.8 Gyr ( Anthony- Twarog, Atwell & Twarog 2005). 
For an age of 3.8 Gyr all NGC 2243 binaries with periods shorter than a few 
days are expected to have circularized orbits (Mathieu 2005). The adjustable 
parameters in the orbital solution were the velocity semi-amplitudes {K\ and 
K 2 ) and the systemic velocity (7). The fit was performed using the GAUSSFIT 
task within IRAF/STSDAS. The derived parameters of the spectroscopic orbit 
are listed in Table 2*. The systemic velocity of NV CMa agrees within the 

tlRAF is distributed by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the NSF. 

^Through this paper we adopt the following values of constants: Rq = 6.95508-B5 km, 
M Q = 1.9891B30 kg, G = 6.67259B - 11 m 3 kg- 1 S - 1 . 
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Table 1: Radial velocities of NV CMa and residuals from the adopted spectroscopic 
orbit 



Phase HJD-2453000 RVi a RV i RV 2 <trv2 (O-C)a (O-C)b 



0.2634 


724.7208 


-67.20 


2 


.77 


192 


.23 


2.69 


-0.81 


0.12 


0.3178 


636.8353 


-56.61 


1 


.50 


182. 


.27 


1.51 


-1.25 


1.40 


0.3375 


636.8587 


-45.65 


1 


.51 


173. 


.65 


2.04 


2.26 


0.36 


0.3557 


636.8803 


-40.76 


1 


.31 


161. 


.51 


1.80 


-1.23 


-3.24 


0.4344 


282.7961 


12.29 


1 


.26 


114. 


.13 


1.18 


2.09 


0.00 


0.5930 


281.7961 


134.08 


1 


.06 


-11 


.81 


1.10 


1.47 


-1.32 


0.6100 


281.8163 


143.47 


1 


.02 


-22 


.17 


1.10 


-0.17 


-0.45 


0.6259 


281.8353 


152.98 


1 


.18 


-31 


.36 


1.14 


-0.13 


0.00 


0.6422 


281.8546 


160.28 


1 


.22 


-40 


.58 


1.26 


-1.59 


-0.30 


0.6582 


281.8736 


170.66 


1 


.31 


-49 


.45 


1.39 


1.21 


-1.45 


0.8310 


633.8797 


173.11 


1 


.26 


-51 


.55 


1.55 


-0.85 


1.03 


0.8401 


633.8905 


170.83 


1 


.35 


-46 


.98 


1.47 


0.63 


1.77 



Table 2: Orbital parameters of NV CMa. 



Parameter 



Value 



P (days) 

To (HJD-244 0000) 

7 (km s _1 ) 
e 

K x (km s- 1 ) 
K 2 (km s" 1 ) 
Derived quantities: 
A sin i (Rq) 
Mi sin 3 i (Mq) 
M 2 sin 3 i (Mq) 
Other quantities: 
<7i (km s _1 ) 
<J2 (km s _1 ) 



1.18851590(fixcd) 
8663.70748(fixed) 
61.70± 0.30 
(fixed) 
128.55±0.52 
130.87±0.55 

6.096±0.018 
1.084±0.010 
1.065±0.010 

1.37 
1.37 



measurement errors with the radial velocity of the cluster as estimated in the 
next subsection. 

2.2 Analysis of Broadening Functions 

We have analysed the spectra of NV CMa using code based on the broadening 
function (BF) formalism (Rucinski 2002). The BF analysis lets us study the 
effects of the spectral line broadening and orbital splitting even for relatively 
complicated profiles of spectral lines. This macroscopic velocity information is 
obtained regardless of the parameters of the photosphere such as the temper- 
ature, pressure etc. It retains strict linearity in reproduction of the individual 
contributions to the broadening profile, i.e. the individual luminosities of com- 
ponents can be simply estimated from the strengths of the respective peaks in 
the broadening profile. 

In the case of detached binary systems the dominating effect on the observed 
broadening of the spectral lines comes from the rotation of the components 
induced by the orbital motion. If we make the approximation of rigid rotation 
of perfectly spherical stars we can integrate the light over the visible hemispheres 
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analytically. The resulting theoretical rotational BF has a shape described by: 

BF rot (v) = A[(l-P)VTlfi + ^P(l-a 2 )]+C (2) 

v - v rad 

a = 

v rot sin I 

where A is a normalization constant, (3 is the linear limb darkening coefficient, 
C is the continuum level, v ra d is the radial velocity of the center of mass of the 
star, v ro t is the linear velocity on the equator of the rotating star, and i is the 
inclination angle between the axis of rotation and the direction to the observer. 

We have extracted BFs from all of the spectra of NV CMa and have con- 
ducted nonlinear least squares fitting of the model profile to the observed BFs 
to measure the radial and projected rotational velocities of the components of 
NV CMa. In both steps of the calculation we have used our programs based on 
procedures from the GNU Scientific Library. 

As it will be shown in Section 3, the components of NV CMa are almost 
spherical and there are no signatures of spot activity. Our model profile is then 
the sum of two theoretical rotational BFs convolved with a Gaussian with a 
standard deviation of 15 km s _1 . The Gaussian represents effects of the instru- 
mental resolution and is a part of the BF method in which modest smoothing 
(adjusted to the width of the spectrograph slit image) is applied to the BFs. 
We used the spectra in the wavelength range from 400 nm to 495 nm. This 
range roughly corresponds to the B passband in the UBV photometric system. 
Since the shape of the rotational BF only weakly depends on the limb darkening 
coefficient, we have adopted a constant /3 = 0.63 as derived from the (B — V, (3) 
relation of van Hamme (1993). The four parameters A, C, v ra d, and v rot smi 
were simultaneously adjusted in the fitting procedure. Figure 2 presents an ex- 
ample of fitting the model to the BF calculated for a spectrum taken at orbital 
phase 0.61. A linear least squares fit to the measured radial velocities gives the 
following orbital solution: 7 = 61.47±0.14 km s" 1 , K x = 130.93 ±0.59 km s" 1 , 
Ki = 128.62 ± 0.42 km s _1 , which lies well within formal errors from the re- 
sult obtained using the cross-correlation analysis. The RMS error of the fit is 
0.93 km s _1 for both components. We measure the projected rotational veloc- 
ities to be Visini = 51.7 ± 1.4 km and V^sini = 52.4 ± 2.5 km s _1 . Note 
that these rotational velocities are very close to those expected if the com- 
ponents are in synchronous rotation in a circular orbit, adopting the period 
from Table 2 and the stellar radii from Table 1 (Vi sini = 52.0 ± 1.3 km s _1 
and V2 sim = 50.1 ± 1.6 km s _1 ). [SMR: Strictly speaking, we should subtract 
quadratically the Gaussian smoothing of sigma — 15 km/s. Then one gets 49.5 
and 50.2 km/s, in almost perfect agreement with the expectations.] 

The integral of a BF profile is proportional to the total flux from a star. Since 
the absolute measurement of this flux requires a perfect match between the tem- 
plate and object spectra we did not attempt to conduct such an analysis. How- 
ever, given that the effective temperatures of the components of NV CMa are 
very similar, the integrals of the BF profiles arc proportional to the monochro- 
matic luminosities of the two components. An analytical integration of the 
profiles leads to the conclusion that the ratio of the luminosities of stars with 
identical spectral types is proportional to the ratio of the products A-v rot smi for 
each of the system components. From a total of 11 spectra with well defined BF 
profiles we obtained Lib/L 2 b = 1.087±0.011 (rms = 0.034), in good agreement 
with the values measured from the photometric observations (see Section 3). 
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Table 3: Radial velocities of stars from the field of NGC 2243 



Name 


HJD-2450000 


RVi 


RV 2 


member? 


H3110 


3281.9127 


59.4(1.0) 




yes 


H4115 


3281.9285 


60.1(1.1) 




yes 


V4 


3635.9006 


117.5(2.4) 


285.9(9.6) 


no 


V5 


3638.8430 


23.50(0.58) 


94.61(0.56) 


yes 


V7 


3638.8217 


96.88(0.80) 


22.8(5.3) 


yes 


V9 


3638.8711 


65.45(0.68) 




yes 



2.3 Velocities of other objects 

Table 2 lists the radial velocities derived for two red giants and four other 
eclipsing binaries from the cluster field. The velocities of H3110 and H4115 are 
close to the velocities of two other cluster red giants observed by Gratton (1982) 
who measured V ra d = +6'2 km ,s _1 and V ra d = +60 km s^ 1 for H4110 and H4209, 
respectively. Based on these four stars we estimate the radial velocity of the 
cluster to be V ra d = +60.4 ± 0.6 km s^ 1 . 

Two velocity peaks are seen in the cross-correlation functions of V4, V5, and 
V7. The mean velocity for the two peaks for V4 is 210.7 km s _1 , well above 
the cluster velocity. We conclude that V4 is a field binary star not related to 
the cluster. In the case of V5, two peaks of very similar shape and height are 
seen in the cross-correlation function. The average velocity of the two compo- 
nents is 59.1 km s _1 . V5 is a very likely member of the cluster with a mass 
ratio close to to unity. Three peaks show up in the cross-correlation function 
obtained for V7, one at 22.8 km s _1 , one at 96.5 km s _1 and a wing to that 
component at 136.8 km s^ 1 . The mean of first two is 59.6 km s~ x , very similar 
to the cluster average velocity. It is possible that the component with highest 
velocity is caused by contaminating light from the close visual companion of 
V7 whose light leaked through the slit during observations (see finding chart 
in KKTS). For V9 only one strong peak was detected in the cross-correlation 
function at a velocity of 65.4 km s -1 , consistent with cluster membership. It 
is possible that this binary was observed close to conjunction, and further ob- 
servations are needed to see if the secondary component can be detected in the 
cross-correlation function. In conclusion, variables V5, V7 and V9 are likely 
members of NGC 2243. These three binaries are promising targets for detailed 
observations aimed at the determination of the parameters of their components. 



3 Light Curve Solution of NV C Ma 

We have analysed the BV light curves of NV CMa obtained by KKTS using 
the Wilson-Deviney model (Wilson & Deviney 1971) as implemented in the 
light-curve analysis program MTNGA§ (Plewa 1988). The mass-ratio of the 
binary was fixed at the spectroscopic value of m^/mi = 0.9825 ± 0.0018. The 
gravity darkening exponents and bolometric albedos were fixed at 0.32 and 0.5, 
respectively. The linear darkening coefficients were adopted for the B and V 
filters from van Hamme (1993) for an assumed metallicity of [Fe/H] = — 0.49 
(Gratton & Contarini 1994; Friel et al. 2002). An interpolation routine in the 
PHOEBE package (Prsa & Zwitter 2005) was used to get values corresponding 
to the adopted effective temperatures of the two components. The light curves 

§MINGA is available at |http:/ /ftp.camk.edu.pl/camk/tomek/Minga/| 
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used in the analysis are shown in Fig. 3. They contain 121 points in the B band 
and 446 points in the V band. Outside of the eclipses, we used normal points 
formed by averaging 3 to 7 individual observations. There is no evidence for 
totality in any of the eclipses. The primary and secondary eclipses have very 
similar depths, implying similar surface brightness and in turn similar effective 
temperatures for the two components. The light curve is symmetric, indicating 
that its shape is not noticeably affected by spot activity. 

The average color index near quadrature is < B — V > m ax= 0.439 ± 0.020. 
The quoted uncertainty includes an external error arising from the photometric 
calibration. Adopting J B( J B-V r )=0.055±0.004 (Anthony-Twarog et al. 2005) we 
obtain an unreddened color index at maximum light of (B — V)o = 0.384. Using 
the empirical calibration of Ramirez & Melendez (2005) we get an effective 
temperature of the primary component of NV CMa of T\ — 6522 ± 129 K. The 
uncertainty includes the formal uncertainty in the temperature calibration as 
well as the uncertainty of the color index. As it is shown below, the effective 
temperatures of the components do not differ by more than about 30 K and the 
color index at maximum light can be safely adopted as the color index of the 
primary component. 

The following parameters were adjustable in the light curve solution: the 
orbital inclination i, the non-dimensional potentials f^i and f2 2 , the effective 
temperature of the secondary T 2 , and the relative luminosity of the primary 
Li(V;B). For a fixed value of the mass ratio q the potentials f^i and f2 2 directly 
determine the relative radii of the components r\ and r 2 . In the following 
discussion we list "equal volume" mean radii of the components. Our finally 
adopted solution (see below) implies that for both components of NV CMa 
the difference between "polar" and "point" radii amounts to about 2% An 
unconstrained light curve solution obtained with MINGA is listed in Table 4. 
It is worth noting at this point that the errors returned by MINGA take into 
account correlations between all fitted parameters. One should note the rather 
large uncertainties of the derived relative radii and luminosities. This is not 
an unexpected result given the partial eclipses (Irwin 1962). An additional 
complication is that the effective temperatures of the components of NV CMa 
are very similar. 

The accuracy of the fit can be significantly improved by using information 
about the light ratio of the two components derived from the spectroscopic data. 
A grid of solutions was calculated with a fixed value of fix (this is equivalent 
to fixing the value of the radius <r\ >). The result is presented graphically in 
Fig. 3 which shows the calculated values of Lib /^2B and < r 2 > as a function 
of the assumed value of <r\ >. There is a strong anti-correlation between the 
calculated value of < r 2 > and the assumed value of < r\ > . The solutions fulfill 
the condition Lib/L2b — 1.087 ±0.011 only for a very narrow range of radii, 
from < ri >= 0.2001 to < n >= 0.2014. However, the uncertainty in < n > 
obtained this way is severely underestimated. To get a realistic estimate of the 
errors in the fitted parameters one has to take into account correlations between 
them. To attain that goal we derived solutions for 3 fixed values of Lib l^2B 
spanning the range 1.087±0.011. The adjustable parameters were i, fii, f2 2 , T 2 
and (L\v I Liv) ■ The final light curve solution along with the uncertainties of all 
fitted parameters is listed in Table 5. Figure 5 shows the residuals corresponding 
to this solution. 



^Thc "polar radii" is a radius toward the stellar pole while "point radii" is the radius 
toward the Lagrangian point LI of the binary orbit. 
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Table 4: An unconstrained light curve solution 
Parameter Value 



i (deg) 


87.08 ± 0.95 




5.959 ± 0.332 


0,2 


6.141 ± 0.392 


Ti (K) 


6522 (fixed) 


T 2 (K) 


6512 ± 123 


(Liv/Z/2v) 


1.110 ± 0.017 


{Lib/ L,2b) 


1.112 ± 0.024 


<ri > 


0.203 ± 0.013 


<ri > 


0.193 ± 0.014 


rms (V) (mag) 


0.008 


rms (B) (mag) 


0.006 



Table 5: Constrained light curve solution with (L 1b /L 2 b) = 1.087±0.011 
Parameter Value 



i (deg) 


87.09 ± 0.66 




6.010 ± 0.136 


0,2 


6.128 ± 0.155 


Ti (K) 


6522 (fixed) 


T 2 (K) 


6506 ± 75 


(Liv / L2v) 


1.088 ± 0.024 


<ri > 


0.200 ± 0.005 


<r 2 > 


0.193 ± 0.006 


rms (B) (mag) 


0.008 


rms (V) (mag) 


0.006 



4 Absolute properties 

The absolute parameters of NV CMa obtained from our spectroscopic and pho- 
tometric analysis are given in Table 6. The errors in the temperatures include 
all sources of uncertainties. The absolute visual magnitudes My were calculated 
using bolometric corrections derived from relations presented by VandenBerg & 
Clem (2003). The observed visual magnitudes derived from the light curve so- 
lutions are Vi = 17.107±0.023 and V 2 = 17.199±0.024, where the uncertainties 
include the errors of the photometric zero point. For the B band we obtain 
Bi = 17.549±0.021 and B 2 = 17.640±0.021. Figure 6 shows the location of the 
individual components of NV CMa on a color-magnitude diagram at the turnoff 
region of NGC 2243. Note that NV CMa lies on the binary sequence in this clus- 
ter (Bonifazi et al. 1990) while the individual components lie on the sequence of 
single stars, further confirmation of the cluster membership of NV CMa. Using 
these values for the observed and absolute visual magnitudes one obtains an 
apparent distance modulus (to— M)v = 13.25±0.10 and (to — M)v = 13.25±0.15 
for the primary and the secondary components of the binary, respectively, with 
an average value of (to — M)y — 13.25 ± 0.08. This value is in good agree- 
ment with recent determinations of (to — M)y = 13.15 ±0.1 obtained from the 
isochrone fitting method by both Anthony- Twarog et al. (2005; they adopted 
[Fe/H] = -0.57) and VandenBerg et al. (2006; they adopted [Fe/H] = -0.61). 
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4.1 The age 

It is possible to estimate ages of the components of the binary by using theo- 
retical age-luminosity relations. It is worth noting that age-luminosity relations 
based on stellar models are unaffected by uncertainties associated with model 
isochrones relating T e ff and L^oi to color index and absolute magnitude in a 
selected band. In particular, they are unaffected by the way in which models 
treat sub-photospheric convection in low mass stars. In Fig. 7 we show age 
versus luminosity relations based on evolutionary tracks recently published by 
VandenBerg et al. (2006). The tracks were derived from a set of isochrones 
calculated with the program vriso distributed with the model grids. Models for 
[Fe/H] = —0.525 (see the following section) and [cn/Fe] =0.3 were used. Horizon- 
tal lines in Fig. 7 mark L±<jl ranges for a given component. The intersections 
of age-luminosity relations with lines marking 1 a limits of the luminosities give 
limits on the age for a given mass. The 1 a age limits are 4.00 <t\ < 5.8 Gyr 
and 4.1 <t 2 < 6.3 Gyr for the primary and secondary components of NV CMa, 
respectively. The large ranges result mainly from the fact that both compo- 
nents are still in a relatively slow phase of their evolution: as can be seen in 
Fig. 6 they are located about 1 mag. below the turn-off point on the cluster 
color-magnitude diagram. One may also note from Fig. 7 that the errors of the 
luminosities and masses contribute an equal 0.6 Gyr to the total uncertainty in 
the estimated ages. 

Figure 8 shows the time dependence of the radius, also a sensitive diagnostic 
for the age, using the same set of evolutionary models. As before, the solid lines 
correspond to evolutionary tracks for the masses of components of NV CMa. 
The measured radii are indicated by the horizontal lines spanning the range 
±1 a. 

^From Fig. 8 we derive 1 a limits on the ages of the components of NV CMa 
of 3.2 < t\ < 4.6 Gyr and 3.2 < t\ < 4.9 Gyr. These limits are consistent with 
these derived from the age-luminosity relations. The age-radius relations sug- 
gest slightly lower ages in comparison with the age-luminosity relation for both 
components. The overlap of the age estimates implies an age for NV CMa of 
approximately 4.1 to 4.6 Gyr. 

4.2 Metallicity 

The age estimates presented in the previous section suffer from a potential 
systematic error arising from the adopted metallicity of NV CMa. There are 
three modern determinations of the cluster metallicity. Gratton & Contar- 
rini (1994) obtained high resolution spectra of two cluster giants and derived 
[Fe/H] = -0.48 ±0.15. Fricl et al. (2002) used medium resolution spectra of 9 
stars to derive [Fe/H] = -0.49 ±0.05. Finally, Anthony-Twarog et al. (2005) 
employed uvbyCaH[3 photometry to obtain [Fe/H] =— 0.57±0.03. The weighted 
average of these three determinations gives [Fe/H] = — 0.547± 0.025. This value 
is very close to [Fe/H] = —0.525 for which we extracted the model relations used 
above. However, the estimated age of NV CMa is very sensitive to the adopted 
metallicity. For example, using models for [Fe/H] = —0.397 and [Fe/H] = —0.606 
we obtain ages t = 5.425 ±0.025 Gyr and t = 3.85 ±0.45 Gyr, respectively, where 
the very small formal errors just indicate the marginal overlap between the age 
ranges obtained from age-luminosity and age-radius relations. In summary, for 
the masses and ages relevant to the present discussion, the metallicity of the 
analysed stars has to be known with a relative accuracy of a few percent to 
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Tabic 6: Absolute parameters of NV CMa 
Parameter Value 



A (R Q ) 


6.104±0.018 


M 1 (M Q ) 


1.089±0.010 


M 2 (M e ) 


1.069±0.010 


Ri (Re) 


1.221±0.031 


R2 (R e ) 


1.178±0.037 


Ti (K) 


6522±129 


T 2 (K) 


6506±149 


Lboli{Lo) 


2.42±0.23 


Lboh(LiQ 


2.23±0.25 


Myi (mag) 


3.86±0.10 


Mv2 (mag) 


3.95±0.15 



allow meaningful comparison of observational data with the models. This is 
illustrated in Fig. 9 which shows age-luminosity and age-radius relations for a 
star with m= 1.089m© and for 3 values of metallicity. The plotted relations are 
based on models by VandenBerg et al. (2006) for [a/Fe] =+0.3. At a given age 
the separation between relations for close, but different, metallicities exceeds the 
uncertainties of the luminosities and radii of components of NV CMa obtained 
in our analysis presented above. 

5 Discussion and summary 

The analysis of photometric and spectroscopic observations of the eclipsing bi- 
nary NV CMA has allowed us to determine the orbital parameters and physical 
properties of the component stars. Our determinations have formal uncertain- 
ties of 1% in the masses and 3% in radii. Uncomfortably large uncertainties in 
the radii result from the degeneracy of the light curve solution for this partially 
eclipsing system. Comparison with model tracks for [Fe/H] = —0.525 give an 
age of 4.1 to 4.6 Gyr for the binary and consequently for the cluster. This de- 
termination suffers from a substantial systematic error related to the uncertain 
metallicity of the cluster. For the relevant range of stellar masses and ages an 
uncertainty in the metallicity of 0.1 dex leads to an uncertainty in the estimated 
age of about 0.8 Gyr. 

There is still room for obtaining a better determination of the parameters of 
NV CMa, including its age. First of all, it is straightforward to derive masses of 
the components with an accuracy of 0.5% or even better. We used 12 spectra 
only and additional radial velocity data obtained near quadrature will lead to 
an improvement in the estimates of the masses of the components. It is also 
possible to obtain better light curves than these used in the present analysis; 
an improvement is possible in both the quality and the phase coverage. In 
particular, our photometry for the B band had poor coverage inside the pri- 
mary eclipse. Improved light curves would better constrain the parameters of 
the components obtained in the light curve solution. And last, but not least, 
the age determination of NV CMa and of other cluster binaries would benefit 
enormously from accurate metallicity determinations, preferably based on high 
resolution spectroscopy of an extended sample of cluster member stars. 

The analysis of single spectra obtained for four other eclipsing binaries in 
the cluster field indicates that three of them are radial velocity members of 
NGC 2243. The fourth star is definitely a non-member. Further observations 
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of member binaries would allow a better constraint on the cluster age as well as 
a test of evolutionary models of low-mass stars with [Fc/H] « —0.5. 




Fig. 1. Spectroscopic observations and adopted orbit for NV CMa. 
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Fig. 3. Phased BV light curves of NV CMa. 
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Fig. 4. Dependence of the luminosity ratio [L\bI^ib) an d relative radius < T2 > on the 
assumed radius < n >. 
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Fig. 5. The residuals for the fit corresponding to the light curve solution listed in Table 5. 
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Fig. 6. Color-magnitude diagram for the turnoff region of NGC 2243. The location of NV CMa 
is marked with an open triangle. The squares denote the locations of individual components 
of the binary. 
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Fig. 7. Theoretical age-luminosity relations for masses Mi = 1.089 ± O.OIOM0 (lower) and 
M2 = 1.069 ± O.OIOMq (upper). Horizontal dashed lines mark 1<t ranges for the observed 
luminosities of the components of NV CMa. The solid lines represent the mass ± 1 a. 
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Fig. 8. Theoretical age-radius relations for masses Mi = 1.089 ± O.OIOMq (lower) and M2 = 
1.069 it O.OIOMq (upper). Horizontal dashed lines mark la ranges for the observed radii of 
the components of NV CMa. The solid lines represent the mass ± 1 a 
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Fig. 9. Theoretical age-luminosity and age-radius relations for a star with mass M = 1.089Mq 
and for mctallicity [Fe/H] = —0.606 (upper curve), [Fc/H] = —0.525 (middle curve) and 
[Fe/H] = -0.397 (lower curve). 



